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iNTRODUCTiOH 


Tills  import  presents  the  results  of  a  stui'iy  of  the  effects  of  terrain  and  thermal  strarifi'  a- 
rioa  upon  the  flew  in  the  surface  boundary  layer  and  the  planetary  boundary  layer  of  the  atmos¬ 
phere.  The  in vestigations  carried  out  were  both  empirical  and  theoretical,  the  observational 
material  consisting  of  the  vrinds  observed  in  the  lowest  5  km  of  the  atmosphere  over  Santa 
Monica,  California.  Here  a  persistent  strong  temperature  inversion  exists  many  days  of  the 
year  between  approximately  500  and  1000  m,  with  superadiabatic  lapse  rates  below  (daytime) 
and  subadiabatic  lapse  rates  above.  The  surrounding  terrain  includes  a  bay  shaped  coast¬ 
line  (partly  mountainous  and  partly  plain),  the  coastal  plains  and  valleys  backed  up  by  moutain 
ranges  2  to  3  km  in  height. 

The  report  consists  of  three  parts.  In  Part  I  the  Santa  Monica  rawin  data  is  analysed  and 
the  features  of  the  mean  flov/  and  the  diurnal  fluctuations  about  it  are  related  to  the  inversion 
and  the  terrain  features.  The  character  of  the  basic  flow  below,  within  and  above  the  inver¬ 
sion  is  described  and  discussed.  The  nature  of  the  diurnal  fluctuations  about  this  mean  flow 
is  scrutinized  in  detail.  It  is  decom;  osed  into  two  oscillations  possessing  different  orienta¬ 
tions  in  the  horizontal,  different  distributions  in  the  vertical  and  different  phases  in  both  the 
time  and  vertical  coordinates.  Both  oscillations  are  associated  with  a  major  terrain  feature 
and  the  influence  of  each  feature  is  seen  to  be  related  to  the  position  of  the  temperature  in- 
verson.  ^ 

Part  II  of  the  report  consists  of  a  theoretical  analysis  of  the  effects  of  an  inversion  layer 
and  a  sea-breeze  on  the  stationary  flow  in  a  planetary  boundary  layer.  It  is  shown  that  an 
intense  temperature  inversion  virtually  acts  as  a  boundary  surface  which  divides  a  planetary 
boundary  layer  into  two  layers.  In  the  lower  layer  the  flow  primarily  is  thermally  driven  where¬ 
as  in  the  upper  layet  the  flow  is  similar  to  that  in  a  boundary  layer  of  a  homogeneous  fluid. 

The  inversion  layer  moves  in  the  direction  of  the  resultant  of  the  stresses  acting  on  the  layer 
by  the  fluids  above  and  below  it.  A  theoretical  model  is  constructed  and  checked  against  the 
observed  flow  at  Santa  Monica. 

Part  Hi  presents  a  theoretical  discussion  of  the  very  lowest  layer  of  the  atmosphere,  the 
surface  boundary  layer,  and  its  response  to  different  thermal  stratifications.  An  analysis  is 
made  of  some  characteristics  of  the  steady  turbulent  transfer  in  the  boundary  layer  of  a  strati¬ 
fied  fluid.  The  effect  of  the  heat  flux  on  the  variation  of  the  mixing  length  and  the  flux 
Richardson  number  with  height  is  determined.  The  velocity  and  temperature  profiles  are  de¬ 
rived.  It  is  found  that  for  a  constant  free  stream  velocity  an  upward  heat  flux  increases  the 
friction  velocity,  whereas  a  downward  heat  flux  decreases  tiie  friction  velocity.  The  lower 
limiting  value  of  the  flux  Richardson  number  is  found  to  be  —0.5,  which  together  with  the  upper 
limiting  value,  0.5,  obtained  by  Townsend,  gives  the  range  of  the  flux  Richardson  number. 
Velocity  profiles  for  the  non-neutral  conditions  converge  in  the  higher  level  towards  the  pro¬ 
file  for  the  neutral  condition,  a  characteristic  which  agrees  with  the  classical  velocity  profiles 
observed  by  Thornthwaite  and  Kaser. 
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PART  I 


WIND  STRUCTURE  IN  THE  LOWEST  5  KM  ABOVE  SANTA  MONICA,  CALIFORNIA 

by  James  G.  Edinger 

1.  Introduction 

This  section  describes  and  analyses  the  summertime  wind  structure  in  the  lowest  5  km  of 
the  atmosphere  above  Santa  Monica,  California  and  attempts  to  determine  the  influence  of 
complex  orography  and  a  persistent  inversion  on  the  flow  across  the  coastline. 

2.  Geography 

The  geometry  of  the  Santa  Monica  area  deviates  in  several  repects  from  a  flat  plain  with 
a  straight  coastline.  As  seen  in  figure  1-1  the  coastline  is  bay  shaped  consisting  of  an  east- 
west  shoreline  at  the  base  of  the  Santa  Monica  Mountains  and  a  north-northwest  south-southeast 
oriented  shoreline  at  the  western  terminus  of  the  Los  Angeles  coastal  plain.  The  dominating 
orographic  feature  inland  is  a  complex  of  mountain  ranges  about  25  miles  north  of  Santa  Monica 
oriented  roughly  east-west  and  rising  on  the  average  to  2  km,  some  peaks  exceeding  3  km. 
Beyond  these  ranges  lie  the  high  deserts  and  the  southern  end  of  the  San  Joaquin  Valley. 


3.  Data 

The  data  used  are  the  checked  radjosonde  winds  aloft  reports  for  Santa  Monica  for  July 
and  August,  1957  and  1958.  Chosen  for  study  were  those  days  on  which  a  strong  inversion 
was  reported  by  the  1600  PST  radiosonde  and  on  which  the  wind  at  the  3  km  level  was  from  the 
40 “sector  betweeyl66“  and  206.°  About  one  out  of  four  days  met  both  requirements  providing 
a  sample  size  of  29  days.  The  subsequent  discussion  is  based  on  resultant  winds  computed 
from  this  data  for  the  surface,  150,  300,  500,  1000,  1500,  2000,  2500,  3000,  4000,  and  5000  m 
levels  for  the  four  observation  times  of  0400,  1000,  1600  and  2200  PST. 


Figure  1-2  gives  the  mean  position  and  strength  of  the  inversion  for  1600  PST  on  these 

days.  The  inversion  layer  in  the  mean  extended 
from  328  m  to  820  m  with  an  increase  of  tempera¬ 
ture  from  the  bottom  to  top  of  8.7 °C.  The  layer 
below  the  inversion  was  superadiabatic  in  the 
mean  with  a  lapse  rate  almost  double  the  dry 
adiabatic  rate.  The  layer  above  was  subadia- 
batic  with  a  lapse  rate  about  two  thirds  of  the 
dry  adiabatic  rate. 
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Figure  1-2 

Mean  temperatiire  .soumling  at  Santa  Monica,  l600  PST 
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Figure  t-3  gives  the  height*  hodograph 


of  the  resultant  winds  obtained  when  data  from 
all  four  observations  times,  0400,  1000,  1600 
and  2200  PST,  are  combined  to  provide  a  repre¬ 
sentation  of  the  basic  flow. 

Below  the  inversion,  in  the  layer  from  the 
surface  up  to  300  m,  the  wind  is  directed  inland 
approximately  normal  to  the  Santa  Monica  coast¬ 
line  at  about  2  m/sec,  the  result  of  a  strong  day¬ 
time  sea-breeze  added  to  a  weak  nighttime  land- 
breeze.  Above  the  inversion  the  wind  velocity 
increases  systematically  from  1.5  m/sec  at  1000 
m  to  7.5  m/sec  at  5000  m  from  the  southerly 
quadrant  and  turns  progressively  to  the  right  with 
increasing  height.  The  southerly  direction  of 
the  flow  at  these  levels  compares  favorably  with 
the  direction  of  the  geostrophic  wind  as  given  on 
the  700  mb  charts  for  these  29  days.  On  most  of 
these  maps  a  north-south  trough  is  indicated 
west  of  the  Southern  California  coast  separating 
the  sub-tropical  high  in  the  Eastern  Pacific 
Ocean  from  a  high  pressure  system  over  the  con¬ 
tinent. 

The  wind  velocity  in  the  inversion  layer, 

300  m  to  1000  m.  has  little  change  in  magnitude 
from  bottom  to  top  but  its  direction  changes  by 
almost  70°  indicating  that  this  extremely  stable 
layer  is  the  region  of  marked  directional  shear 
between  lower  west-southwesterly  and  upper 
southerly  wind  regimes. 


Height  hod'graph  averaged  for  all  observation  times, 

Santa  Monica 


5.  Diurnal  fluctuations 


Figure  1-4  presents  the  diurnal  variations  of  the  velocity  about  this  basic  flow  in  the 
form  of  time  hodographs  for  each  level  from  the  surface  up  to  the  5  km  level.  Each  hodograph 
consists  of  the  four  vectors  representing  the  velocity  relative  to  the  basic  flow  at  0400,  1000, 
1600  and  2200  PST.  These  four  vectors  are  plotted  at  the  end  point  of  the  velocity  vector  of 
the  basic  flow  so  that  the  relative  orientation  and  magnitude  of  the  fluctuation  and  the  basic 
flow  are  apparent. 


♦  Two  types  of  hodographs  will  be  discussed:  (a)  one  showing  the  variation  of  velocity  with 
height  and  (b)  one  showing  the  variation  of  velocity  with  time  of  day.  To  distinguish  the 
two  the  modifiers  "height"  or  "time"  will  be  used. 
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1  i'C  inversion  inycr  is  seen  to  be  the  seat  of  marked  changes  not  only  in  the  basic  flow 
I'u!  i.i  i;w;  cliarac'iT  of  '.he  diurnal  oscillation  as  well.  Whereas  the  mean  flow,  as  mentioned 
above,  changes  direction  abruptly  here,  the  magnitude  of  the  diurnal  oscillation  suffers  an 
abrupt  climunition  in  magnitude. 

In  the  lowest  500  m  the  diurnal  oscillation  of  wind  velocity  is  largest,  about  5  m/sec  in 
magnitude  along  the  60°— 240'’  axis  with  maximum  on-shore  (240°)  component  at  about  1600 
PST.  Above  500  m  the  magnitude  of  the  60°— 240°  oscillation  decreases  rapidly  reaching  zero 
at  1500  m.  Above  1500  m  its  phase  is  reversed  (maximum  on-shore  component  at  0400  PST) 
and  its  manitude  increases  to  a  new  maximum  of  2  m/sec  at  2500  m.  Above  2500  m  the  magni¬ 
tude  of  the  oscillation  decreases  to  about  0.5  m/sec  at  4000  and  5000  m,  the  phase  remaining 
unchanged. 

This  oscillation  in  the  60°— 240°  vertical  plane  has  the  characteristics  generally  ascribed 
to  a  sea-land  breeze  regime.  However,  the  observed  fluctuations  in  velocity  are  not  confined 
entirely  to  the  60°— 240°  vertical  plane.  Components  of  motion  outside  this  plane  become  more 
pronounced  as  height  increases.  At  the  3,  4  and  5  km  levels  the  maximum  oscillation  is  more 
nearly  along  the  180°— 360°  axis  than  the  60°— 240°.  Furthermore,  the  oscillation  in  this 
north-south  plane  apparently  has  a  phase  that  lags  about  6  hours  behind  the  60 °— 240°  oscil¬ 
lation.  Figure  1-5  contrasts  the  observed  changes  in  velocity  from  0400  to  1600  PST  with 
those  from  1000  to  2200  PST  at  all  levels.  Whereas  the  velocity  variations  with  the  0400  to 
1600  PST  phase  are  seen  to  be  confined  mostly  to  the  60°-240°  direction,  those  with  the 
1000  to  2200  PST  phase  appear  to  be  more  nearly  oriented  along  a  180° -360°  axis. 

The  observed  velocity  fluctuations,  then,  are  better  described  as  a  superposition  of  two 
oscillations  than  by  only  one.  Figure  1-6  provides  the  observed  magnitudes  of  the  60°— 240° 
and  180°— 360°  oscillations  as  a  function  of  height.  There  are  several  obvious  differences 
between  the  two  distributions:  (1)  The  strength  of  the  lower  branch  of  the  60°— 240 ''oscilla¬ 
tion  is  much  greater  than  that  of  the  lower  branch  of  the  180° -360°  oscillation,  the  upper 
branches  of  both  being  of  about  the  same  strength,  and  (2)  The  level  separting  the  upper  and 
lower  branches  of  the  two  oscillations  are  displaced  500  m  with  respect  to  each  other,  the 
180°— 360°  oscillation  reversing  at  about  2000  m,  the  60°— 240°  at  1500  m. 

If  the  two  oscillations  were  in  phase  timewise,  the  resultant  hodographs  would  be  simple 
lines,  enclosing  no  area.  Since,  however,  there  is  a  time  lag  between  the  two,  the  hodographs 
open  out  into  closed  curves.  Furthermore,  the  displacement  in  the  vertical  of  the  two  oscil¬ 
latory  pattern  with  respect  to  each  other  determines  the  sense  of  rotation  of  •'he  vectors  in  the 
time  hodographs.  If  both  patterns  reversed  their  direction  of  flow  at  the  same  height,  the  sense 
of  rotation  would  be  the  same  at  all  levels.  Since,  however,  the  reversal  takes  place  at  1500 
m  for  the  60°— 240°  oscillation  and  at  about  2000  m  for  the  180°— 360°,  the  sense  of  rotation 
is  reversed  betweeo  thye  two  levels. 

Figure  hi  ehewi  whnr  the  time  hodograph  would  be  for  eueh  level  if  the  veeror  sum  of 
these  two  osellldtlafts,  6  haufs  aut  af  phase  timewise  and  500  m  aut  af  phase  hclghtwlse, 
completely  described  the  motion.  For  comparison  purposes  it  also  displays  the  corresponding 
observed  time  hodographs. 

There  is  general  agreement  between  the  synthesized  and  observed  time  hodographs.  Even 
the  sense  of  rotation  of  the  vectors  is  correctly  described  at  almost  all  levels,  the  predominant 
counterclockwise  rotation  giving  way  to  clockwise  only  in  the  vicinity  of  the  1500-2000  m 


MAGNITUDE  OF  OSCILLATION  in  in/5 

Figure  1-6  Magnitudes  of  the  60°— 240°  and  180°— 360°  oscillations  as  a 

function  of  height 


layer  where  a  180°  shift  in  the  relative  phase  of  the  two  oscillations  results  from  their  500  m 
relative  displacement  in  the  vertical. 

To  reveal  more  completely  the  details  of  the  motion  in  the  60°— 240°  and  180°— 360° 
vertical  planes,  time  sections  for  both  are  presented  in  figures  1-8  and  1-9.  The  obvious  tilt 
of  the  systems  toward  increasing  time  with  increasing  height  indicates  that  the  phase  of 
neither  oscillation  is  actually  constant  with  height.  Oscillations  at  higher  levels  consistently 
lag  behind  those  below.  Similar  behavior  was  observed  in  the  Batavia  sea-breeze  by  van 
Bemmelen,(l)  (1922),  see  figure  1-10.  The  upper  branch  of  the  circulation  over  Batavia  reaches 
its  maximum  speed  about  1.5  hours  after  the  maximum  occurs  in  the  lower  branch.  For  Santa 
Monica  this  lag  is  2.5  hours  for  the  60°— 240°  circulation  and  3.5  hours  for  the  180°— 360° 
circulation.  The  lag  between  the  60°— 240°  and  180°— 360°  circulations  at  Santa  Monica  is 
seen  ti  he  closer  to  3  hours  than  to  6  hours  as  suggested  by  the  simplified  analysis  above. 

The  general  dimensions  of  the  diurnal  circulation  above  Batavia  and  Santa  Monica  are 
about  the  spme.  The  maximum  return  flow  aloft  occurs  at  2  km  in  the  former  and  at  3  km  in 
th^  latter.  Although  Sknta  Monica’s  circulation  is  deeper  it  is  only  about  half  as  strong. 
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Figure  1-10 

Time  section  of  Batavia  sea-breeze  in  vertical 
plane  normal  to  coast,  (speed  in  m/sec) 


Batavia's  upper  and  lower  Hows  have  niaxi- 
raum  speeds  of  3  and  7  m/sec  respectively 
concra.sted  with  1  and  3  ni/sec  for  Santa 
'Monica.  At  both  cities  the  maximum  depth 
of  the  circulation  is  reached  in  late  afternoon. 
As  interesting  point  of  difference  is  revealed 
by  considering  the  heights  of  occurence  of 
the  maximum  sea-breeze:  at  Santa  Monica  it 
is  found  at  the  surface  and  at  Batavia  at  about 
150  m  above  the  surface.  However,  it  is 
likely  that  the  speed  of  Santa  Monica’s  sur¬ 
face  sea-breeze  is  exaggerated  because  the 
observation  is  taken  at  the  top  of  a  bluff  that 
faces  directly  into  the  sea-breeze.  Pilot 
balloon  observations  taken  at  U.C.L.A.,  a 
few  miles  inland,  indicate  a  maximum  well 
above  the  surface  as  shown  in  figure  1-11, 
the  data  indicating  the  flow  beneath  the  in¬ 
version  to  be  a  modified  channel  flow  with  a 
maximum  speed  a  little  below  the  middle  of 
the  layer,  not  at  the  ground. 


WIND  SPEED  in  m/t 
Figure  1-11 

Speed  of  sea-breeze  versus  height  in  marine  layer, 
U.C.L.A. 


6.  interpretation 

The  foregoing  description  of  the  flow 
across  the  coastline  at  Santa  Monica  has 
several  unique  features  that  distinguish  it 
from  observed  patterns  at  other  coastline  loca¬ 
tions  as  well  as  from  patterns  prescribed  by 
the  various  sea-breeze  theories.  The  attempt 
is  made  here  to  relate  these  distinguishing 
features  to  geographic  and  atmospheric  con¬ 
figurations  peculiar  to  the  Santa  Monica  area. 

The  mean  flow  (diurnal  variations  aver¬ 
aged  out)  has  a  unique  structure  apparently 
related  to  the  persistent  and  strong  inversion 
above  Santa  Monica.  Below  the  inversion  a 
monsoon  flow  exists,  on-shore  parrallel  to 
the  Santa  Monica  mountains.  Its  profile  re¬ 
sembles  channel  flow  but  its  direction  shifts 
gradually  toward  the  direction  of  the  flow 
above  it.  In  the  inversion  layer  the  direction 


of  flow  changes  rapidly  resulting  in  a  strong  directional  sheur. 


Above  the  inversion  the  flow 


turns  slowly  to  the  right  in  a  manner  approximating  an  Ekman  spiral,  finally  assuming  the  direc¬ 
tion  of  the  geosrophic  wind  at  heights  above  2.5  km.  In  the  next  section  of  this  report  Kao 
(1959)  develops  a  theoretical  model  for  steady  flow  in  a  layer  of  air  bisected  horizontally  by  a 
strong  temperature  inversion.  His  model  reproduces  the  major  features  of  this  observed  mean 
flow  above  Santa  Monica. 


C  .  ;i;''  Jiuii!;’-  v:\ri.'‘Uoii;j  in  ri;»-  flov.’  pattern  no  theoretical  description  has  I'.een 

att'-mf'ted  a.iihiM,;;';  the  nat'.ire  oi  certain  re!a; ion-ship;;  between  terrain  features  and  the  observed 
atmospiieric  osciilations  are  suggested  by  cite  data  rxnd  may  provide  fruitful  topics  for  later 
theoretical  investigations.  The  fact  that  the  somewhat  complex  observed  oscillatory  pattern 
can  be  approximated  by  two  superposed  simple  o.scillarions  suggests  that  in  addition  to  the 
area  of  strong  differential  heating  along  the  coastline  there  must  be  somewhere  else  another 
area  almost  as  effective  in  producing  a  diurnal  oscillation  of  flow. 


The  oscillation  along  the  60 240°  axis  has  the  characteristics  of  the  conventional  sea- 
land  breeze  at  a  coastline.  However,  the  coastline  is  strongly  curved  in  the  vicinity  of  Sapta 
.Monica.  Consequently,  it  is  not  clear  just  what  orientation  the  sea-land  breeze  regime  should 
take.  The  fact  that  it  moves  almost  perpendicular  to  the  shoreline  in  the  immediate  vicinity  of 
Santa  Monica  may  in  large  part  be  the  result  of  the  roughly  east-west  steering  action  of  the 
Santa  Monica  Mountains  several  miles  to  the  north.  The  relatively  dense  500  m  thick  layer  of 
cool  marine  air  lying  beneath  the  strong  persistent  inversion  provides  a  very  stable  configura¬ 
tion.  The  Santa  Monica  mountains,  somewhat  higher  than  500  m  on  the  average,  therefore,  con¬ 
stitute  an  effective  constraint  on  the  direction  of  flow  of  this  layer. 


The  oscillation  along  the  180°— 360°  axis,  on  the  other  hand,  lags  3  hours  behind  the 
60 240°  oscillation  might  be  the  result  of  differential  heating  along  an  east-west  line  at 
some  distance  from  Santa  Monica  and  at  a  height  above  the  steering  influences  of  the  Santa 
Monica  mountains.  Such  heating  is  likely  along  a  different  sort  of  "coastline"  some  20  miles 
north  of  Santa  Monica,  the  line  of  demarkation  that  lies  between  the  top  of  the  marine  layer  and 
the  mountain  ranges  that  rise  above  it.  The  shaded  area  in  figure  1-1  represents  the  land  that 
stands  above  this  "ocean”  of  marine  air,  i.e.  the  terrain  above  the  1500  ft  contour.  It  has, 
roughly  speaking,  an  east-west  orientatior .  Of  course,  even  in  the  absence  of  the  marine  layer 
and  the  strong  inversion  differential  heating  would  take  place  here  between  the  air  in  contact 
with  the  mountains  and  that  at  the  same  level  but  above  the  coastal  plain.  The  inversion,  how¬ 
ever,  by  preventing  the  convective  transfer  of  heat  to  the  air  layers  above  the  marine  layer  in¬ 
tensifies  and  localizes  the  differential  heating  along  a  zone  at  about  500  m  elevation  along  the 
southern  slopes  of  the  mountain  ranges.  Interesting  corroborative  evidence  may  be  the  observed 
fact  that  the  180°— 360°  oscillation  is  displaced  vertically  about  500  m  with  respect  to  the 
60°— 240°  oscillation. 


If  thi.s  upper  "coastline,”  20  miles  inland  of  the  regular  coastline,  is  the  seat  of  the  180°  — 
360°  oscillation,  the  effect  of  the  differential  heating  must  be  propogated  southward  at  a  rate 
of  3.5  m/sec  relative  to  the  ground,  6.5  m/sec  relative  to  the  air  (3  m/sec  headwind)  in  order 
to  arrive  over  Santa  Monica  with  the  observed  phase  lag  of  3  hours.  This  rate  of  propogation  is 
similar  to  that  obtained  from  the  Batavia  data.  There  a  mountain  range  parallels  the  coast  some 
50  to  75  km  inland.  The  phase  lag  at  the  500  m  level  (about  50  km  from  the  coast)  is  about 
1.5  hours.  Assuming  no  mean  wind,  the  monsoon  roughly  parallel  to  the  mountains,  the  rate  of 
propogation  would  be  10  m/sec.  It  appears  that  at  both  Batavia  and  Santa  Monica  the  phase 
lag  with  height  could  be  ascribed  to  effects  originating  on  the  distant  mountain  slopes. 

Just  as  the  Santa  Monica  mountains  steer  the  low  level  fluctuations,  the  inland  mountain 
ranges  apparently  block  them.  This  is  evident  in  the  180°— 360°  component  of  the  fluctuations. 
It  was  pointed  out  in  the  discussion  of  figure  1-6  that  the  lower  branch  of  the  180°— 360°  oscil¬ 
lation  was  much  weaker  than  the  upper  branch.  It  is  significant  that  the  depth  of  the  lower 
branch  just  about  matches  the  average  height  of  the  mountain  range,  approximately  2  km.  Fluc¬ 
tuations  at  the  mountain  slopes  are  constrained  to  move  up  and  down  the  slopes  producing  a 
vertical  fluctuation.  This  vertical  fluctuation  must  have  maximum  amplitude  near  the  ridge  top 
and  decrease  with  increasing  height  above,  producing  a  fluctuating  convergence  of  the  vertical 
component  there.  Such  a  fluctuating  convergence  pattern  would  be  reflected  in  an  increase  in 
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thi;  n!  til'-  iinrizoiUiil  velocity  fiiictuatioiis  in  tlic  upper  flow,  (at  levels  higher  than 

tiie  r'.ijg.cs).  'Ihe  fact  that  the  60'  — 240"'  oscillation  at  heights  above  2  km  is  just  as  large  as 
t!ic  180  '~3()0  ’  oscillation  yet  receives  no  amplification  due  to  mountain  slopes  is  one  of  the 
.several  iiulication.s  rhar  this  oscillation  is  the  more  vigorous  of  the  two.  It  .seems  likely  that 
the  proximity  of  Santa  Monica  to  the  coa.-;:  and  remoteness  from  the  inland  mountain  ranges 
would  prescribe  that  the  60°— 2''i0°  oscillation  be  the  dominant  component. 

One  component  prescribed  by  the  sea-breeze  theories  undetectable  in  the  observations  is 
that  due  to  the  Coriolis  acceleration.  Theory  based  on  flat  terrain  and  a  long  straight  coast¬ 
line  requires  the  time  hodographs  to  have  a  clockwise  rotation.  At  Santa  Monica  they  are 
observed  to  have  the  opposite  rotation  almost  exclusively,  the  major  exception  being  in  the 
layer  between  1500  and  2000  m  where  .he  effect  of  the  lower  branch  of  the  fluctuation  set  up 
by  the  inland  mountain  ranges  combines  with  the  upper  branch  of  that  originating  at  the  coast 
to  produce  a  clockwise  rotation.  It  is  apparent  that  the  influence  of  the  inland  mountain  ranges, 
like  that  of  the  coastline,  is  a  first  order  effect  in  determining  the  nature  of  the  diurnal  flow 
patterns  and  that  the  coriolis  effect  must  be  regarded  as  of  second  order. 


7.  Conclusions 

The  mean  flow  of  air  over  Santa  Monica  in  summer  in  the  presence  of  a  southerly  geostro- 
phic  wind  a:  700  mb  can  be  convenieni'y  subdivided  into  three  layers:  (1)  below  the  inversion 
—  monsoon,  (2)  in  the  inversion  —  strong  directional  shear  and  (3)  above  the  inversion  —  Ekmon 
spiral.  Superposed  on  this  somewhat  complex  mean  flow  is  a  diurnal  sea-land  breeze  oscilla¬ 
tion  consisting  of  a  lower  branch  about  as  deep  as  the  ini  cq  noui.  '••.In  ranges  are  high,  (2000 
m)  with  a  return  flo'^ aloft  of  similar  d  ,'pth.  This  diurna!  ^  .-v-'IIUtion  can  Le  represented  as  the 
sum  of  two  different  oscillations  and  therefore  appears  to  be  the  itsult  of  differential  heating 
effects  in  two  areas;  (1)  along  the  coastline  of  Santa  Monica  Bay  and  (2)  along  the  "shore” 
of  the  "ocean”  of  marine  air,  at  500  m  elevation  along  the  southern  slopes  of  the  inland  moun¬ 
tain  ranges  The  difference  in  orientation,  height  add  proximity  of  these  two  "shorelines”  with 
respect  to  Santa  Monica  apparently  prescribe  the  major  features  of  the  complex  configuration  of 
the  sea-land  breeze  circulation. 

The  inversion  influences  the  flow  patterns  by  virtually  prohibiting  the  vertical  transfer  of 
heat  and  momentum.  It  plays  the  role  of  the  "ground”  for  the  flow  above  500  m  above  Santa  Monica. 
This  mean  flow  which  resembles  the  Ekman  spiral  represents  an  elevated  planetary  boundary 
layer.  The  inversion  plays  a  similar  role  in  the  thermal  circulation  set  up  at  the  mountain 
slopes,  in  this  case,  being  analagous  to  the  sea  surface,  the  mountain  slopes  near  the  500  m 
level  being  taken  as  the  beach.  In  the  case  of  the  thermal  circulation  set  up  at  the  real  coast¬ 
line,  however,  it  forms  neither  top  nor  bottom  of  the  circulation.  Here  it  constitutes  a  kind  of 
in.sulating  layer  interior  to  the  lower  branch  of  the  circulation  where  the  magnitude  of  the 
oscillation  begins  to  decrease  rapidly  with  height. 

The  effect  of  the  coriolis  accelaration  on  the  flow'  pattern  is  evidently  overwhelmed  by  the 
much  larger  influences  of  the  terrain  features.  The  clockwise  turning  of  the  velocity  vector  in 
the  time  hodographs  occurs  only  where  terrain  effects  dictate  this  behavior. 

These  results,  though  suggestive,  are  tentative.  They  depend  for  verification  or  repudi¬ 
ation  upon:  (1)  similar  studies  of  observational  material  from  areas  of  diverse  orography  and 
thermal  stratification  and  (2)  theoretical  investigations  of  disturbances  in  the  field  of  motion 
brought  about  by  differential  heating  on  a  diurnal  schedule  associated  with  terrain  features 
and  coastlines. 
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PART  !! 

STATIONARY  FLOV/  IN  THE  PLANETARY  BOUNDARY 
LAYER  WITH  AN  INVERSION  LAYER  AND  A 
SEA  BREEZE 
S.-K.  Kao 


1.  Introduction 

In  recent  years,  synoptic  analysis  of  the  lower  troposphere  over  the  coastal  plain  of  Cali¬ 
fornia  shows  that  during  the  summer  seasons  a  temperature  inversion  occurs  in  the  layer  between 
500  and  1000  m  on  the  average.  (Neiburger,  Beer  and  Leopold,  1944).  Under  this  inversion 
layer  a  land-sea  breeze  prevails.  On  account  of  the  presence  of  this  inversion  layer,  the  flow 
in  the  planetary  boundary  layer  deviates  greatly  from  that  in  a  homogeneous  fluid.  A  detailed 
description  and  discussion  of  the  diurnal  wind  distribution  at  various  levels  over  Los  Angeles 
Basin  has  been  given  by  Edinger  (1959).  It  is  the  objective  of  this  note  to  examine  the  effect 
of  a  temperature  inversion  and  a  sea-breeze  on  the  stationary  flow  in  the  planetary  boundary 
layer  as  a  first  step  to  the  understanding  of  the  mechanism  of  this  system. 


2.  Quosi-stationary  wind  distribution  in  the  planetary  boundary  layer 


The  presence  of  a  sea-breeze  in  the  lower  layer  of  the  planetary  boundary  layer  makes  it 
difficult  to  attain  a  stationary  flow.  As  a  first  appro-ximation  to  a  stationary  state  and  to  bring 
out  the  effect  of  the  sea-breeze,  we  use  the  observation  taken  at  the  time  when  the  sea-breeze 
approaches  its  maximum.  To  bring  out  the  effect  of  the  wind  at  the  geostrophic  level  we  con¬ 
sidered  the  winds  at  this  level  occurring  in  the  40*’  sector  from  l66°to  206°  degrees.  In  so 
doing  we  obtain  a  sample  of  data  of  29  cases  from  the  rawin  soundings  taken  during  July  and 
August,  1957  and  1958,  at  Santa  Monica  Airport,  California.  The  average  temperature,  Fig.  1-2 
shows  an  intense  inversion  layer  occurring  at  a  height  between  328  m  and  820  m.  The  planetary 
boundary  layer  consists,  therefore,  of  three  layers,  the  inversion  layer,  and  the  layers  below  and 
above  it.  The  distribution  of  the  average  winds  at  various  levels  is  represented  by  large  dots 
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in  Fig.  2-1.  Two  characteristic  features  appear 
in  this  hodograph;  (1)  In  the  lower  layer,  say 
below  1  Km,  the  wind  turns  toward  the  left  with 
height,  whereas  in  the  upper  layer  it  turns  to¬ 
ward  the  right.  (2)  In  the  lower  layer  the  speed 
of  the  wind  first  increases  with  height,  then  de¬ 
creases  to  a  minimum  at  a  height  about  1  Km; 
in  the  upper  layer  wind  speed  generally  increases 
with  height. 

It  is  interesting  to  note  that  in  the  lower 
layer  the  hodograph  bears  no  resemblance  to 
the  classical  wind  distribution  in  a  planetary 
boundary  layer  of  a  homogeneous  fluid  (Ekmaii; 
1902;  Akerblom,  1908;  Taylor,  1916;  Rossby, 

1932;  Rossby  and  Montgomery,  1935).  The 
turning  of  the  wind  toward  the  left  instead  of 
right  with  height  indicates  that  the  Coriolis 
force  is  definitely  small  as  compared  with  other 
forces.  Furthermore,  analysis  of  winds  for 
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variou.s  case;-;  ,slio\v.s  th;u  strong  velocity  shears  exist  between  the  inversion  layer  and  its 
neighboring  fluid.  It  is  clear  that  the  occurrence  of  these  shearing  stresses  is  due  to  the 
intense  stability  of  the  inversion  layer  which  virtually  acts  as  a  boundary  and  prevents  the 
transports  of  both  momentum  and  heat  from  one  side  to  the  other.  The  flow  in  the  lower  layer 
therefore,  is  primarily  governed  by  the  sea-breeze  subjected  to  the  boundary  effects  of  the  in¬ 
version  layer  and  at  the  surface  of  the  earth.  It  is  known  both  from  observations  and  theories 
(Van  Bemmelem,  1922;  Haurwitz,  1947;  Schmidt,  1947;  Defant,  1950;  Pierson,  1950;  Pearce, 
1955;  Smith,  1955)  that  above  a  sea-breeze  a  comparatively  weak  land-breeze  occurs.  On  ac¬ 
count  of  the  presence  of  an  inversion  layer  the  effect  of  land-breeze  in  the  upper  layer  would 
be  greatly  reduced  and  that  the  flow  in  the  upper  layer  is  similar  to  that  in  the  boundary  layer 
of  a  homogeneous  fluid. 


3.  Theoretical  Mode! 

To  construct  a  simple  model  bringing  out  the  major  effects  of  an  inversion  layer  and  a 
sea-breeze  on  the  flow  in  the  planetary  boundary  layer,  we  idealize  the  inversion  layer  to  a 
surface  of  discontinuity  in  temperature.  On  account  of  its  extreme  stability,  the  surface  of 
discontinuity  may  be  considered  as  a  boundary  surface  which  isolates  the  transfers  of  momen¬ 
tum  and  sensible  heat  from  one  side  to  the  other  of  the  surface.  The  planetary  boundary  layer 
consists,  therefore,  of  two  layers,  one  above  and  the  other  below  the  surface  of  discontinuity. 
In  the  upper  layer  the  flow  is  controlled  by  the  pressure,  Coriolis  and  frictional  forces,  where¬ 
as  in  the  lower  layer  we  assume,  fo.  simplicity,  that  the  flow  is  primarily  thermally  driven  and 
that  the  Coriolis  effect  may  be  neglected.  The  equations  of  motion  for  the  upper  and  lower 
layer  are  respectively 

Ki^— (ui +ivi)  -  if  (u^ +ivi;  =  .L(^^  +  i^)= -if(u  +v  )  (1) 

dz2  Pi  Bx  By 

K2^(u2MV2)-  JL(!L2.i!P2)=o 

dz^  p2  Bx  By 


where  K  is  the  coefficient  of  eddy  viscosity,  p  the  density,  p  the  pressure,  f  :=2  fl  sin<^> 
is  the  Coriolis  parameter;  subscripts  1  and  2  are  respectively  referred  to  the  upper  and 
lower  layer,  the  x,y,z  axes  are  respectively  oriented  toward  the  east,  north,  zenith.  V  =u+iv 
and  Vg  =  Ug  +  iVg  are  respectively  the  horizontal  and  geostrophic  velocities. 

Let  the  horizontal  axes  be  coincided  with  the  surface  of  discontinuity  which  is  assumed 
to  be  horizontally  oriented.  The  kinematic  boundary  condition  requires  the  velocity  to  be  con¬ 
tinuous  at  the  surface  of  discontinuity  of  which  the  velocity  Is  denoted  by  Vd  «=  ud+ivd.  V?e 
have 


at  z  =  0  :  ui  =  U2  =  uj  ,  vi=  V2  =  vj  (3) 

Furthermore,  we  assume  that  the  surface  of  discontinuity  acted  on  by  the  air  from  above  and 
below  moves  in  the  direction  of  the  resultant  of  the  velocity  shears  at  this  surface 


Bui  Bu-1  Bvi  Bv 

at  z  =  o:  uj  =  Cj,  -  2^,  vj  =  — 

Bz  Bz  Bz  Bz 


2., 


(4) 


where  c^  is  a  constant. 
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I'liL-  Other  boiuKlnt'.'  conditions  arc 


at  00  :  u  =  u  ,  V  =  v 

o  o 


at  2  =  -H  :  u  =  V  =  0 


where  H  is  the  height  of  the  surface  of  discontinuity. 


The  general  solution  of  (1)  is 


-{I  *i)  1,2.  (1  +i)  );z 

ui  +  ivj  .  Ug  ..Vg  +  cje  .  C2  e  ^2) 

where  ^  =  ({/2K-^)''^  ,  cj^  and  C2  are  two  constants. 

Applying  the  boundary  condition  (5)  to  (7),  we  have  C2  =0.  Eq.  (7)  becomes 

ui +ivi  =  Ug +iVg  +  cje'(^  (8) 

For  the  convenience  of  later  discussion  let  H  be  the  height  at  which  the  wind  blows  parallel  to 
the  isobars.  We  have  4  =  tt/H  and,  therefore  ffH/w)^. 

If  we  now  apply  the  kinematic  boundary  condition  (3)  to  (8),  we  have 

Ci  =  {(uj-ug)  M(vj-vg)  I  (9) 

Substitution  of  (9)  in  (8)  gives 


for  2  >  0 


For  simplicity,  we  assume  that  the  horizontal  pressure  gradient  is  constant  in  the  lower 
layer.  It  can  be  shown  that  the  solution  of  (2),  which  satisfies  the  boundary  conditions  (3) 
and  (6),  takes  the  form 

for  0  >  z  >  -H. 


Let  us  define  the  sea-breeze  velocity,  =  Ug  +  ivg,  as  the  velocity  in  the  lower  layer 
if  the  surface  of  discontinuity  is  at  rest.  Putting  u^j  =  vj  =0  in  (11)  we  have 


y^Kjdx  3y  I'  2  '  '2  M. 


2/^2^ 
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(13) 


F'.jrth.t-r'nore,  let  the  'iiaxirnuni  .eca-breeze  velocity,  at  z  -  -  H'/7.  ,  be 


=  ‘‘s  ^  '^s  = 


(1P2  .  1  ;iP2.) 


8  K2  'lx  3y 

The  pressure  force  in  the  lower  layer  may,  therefore,  be  expressed  in  terms  of  the  maximum 

sea-breeze  velocity,  Vg  .  Elimination  of  the  pressure  gradient  from  (11)  and  (13)  gives 

^  2 

U2  +iv2  «  H  -  (1  +A  z)  1  (Ug  +ivs)  +  .+— )  (ud  +ivd) 

I  '  H  '  J  '  H/  (14) 


The  velocity  distribution  of  the  fluid  below 
the  surface  of  discontinuity  in  relation  to  the 
maximum  sea-breeze  velocity,  Vg  ,  and  the 
velocity  of  the  surface  of  discontinuity  may 
be  seen  in  Fig.  2-2.  The  quantity  remaining 
to  be  determined  in  (10)  and  14)  is  the  velocity 
of  the  surface  of  discontinuity.  This  can  be 
done  using  the  condition  (4).  We  obtain  from 
(10)  and  (14)  the  velocity  shear  at  the  surface 
of  discontinuity. 


Schematic  representation  of  sea-breeze  profile 
Figure  2-2  and  velocity  of  the  surface  of  discontinuity 


|4ui+ivi)| 


=  -<1 +i).^{(Ujj  -  Ug)  +  i(vj-v^) } 


z  =0 


H 


_(u2+iv2)  = 

z  =  o 

Substitution  of  (15)  in  (4)  gives 


-i{(udMvj)  -4(Es  ^ivg)} 
H 


(15) 


(Ujj+ivj)  =  Cj,-T{(l +i)  (u  +iv  )  +  l.(Ug +iv5) -(1 +i. +i)  (uj +ivjj)} 

H  -V  TT 

Solving  for  (ud  +  ivd),  we  obtain 

((2+ii_  )  +  i  H  1  (ug+iv  )  +i-|(l +iL) -iUiTs +ivs) 

^  ^  J  _ 

1  +  A  +i_+JL.)2 

\  -FT  C„Tr  ./ 


“d  +  i^d 


(16) 


It  can  be  shown  that  c^  is  of  the  order  of  magnitude  of  H/n  .  To  estimate  the  velocity  of 
the  surface  of  discontinuity  we  put  c^  =  H/tt  in  (I6). 


UJ+IVJ  . 


i(l  +27r)  +i  }  (Ug  +  ivg)  +  4  {(1  +-L  ^  ^  ^  (“s  *  ’^s) 


TT  [  1  +(2  ] 


(17) 
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Subsutution  of  (17)  in  (10)  and  (14)  gives  respectively  the  velocity  distribution  in  the  upper 
and  lower  layers 


t  1 

u j  r  rvj  =  ,  1  +  A 

77 

v 


1  r  ~( 3  +  3  77  +  -L.)  +  i 


1  ^(2-i-)' 


"(1  +i)C2 


1  <-8  *  ‘ 


For  a  given  pair  of  geostrophic  wind  and  maximum  sea-breeze  velocities  the  velocity  of 
the  inversion  layer  may  be  predicted  by  (17),  and  the  velocity  distribution  in  the  planetary 
boundary  layer  may  be  calculated  from  (18)  and  (19). 

The  solid  curve  in  Fig.  2-1  is  computed  from  (18)  and  (19)  in  which  the  geostrophic  velo¬ 


city  is  estimated  from  the  wind  at  2500  m  whereas  Ug  +ivg  is  estimated  by 

(u+iv)2=. 

=  2'^  [(|_  +  ±)  -i]  ^  {77  [  1  +(2  +1)^]  (u  +iv)2.  =  H  [(1  +277)  +i]  (ug  +ivg)  } 

as  shown  in  Fig  2-3,  and  the  wind  at  the  middle  of  the  lower  layer  is  approximated  by  the  wind 
at  500  m. 

4.  Discussion 


It  is  seen  from  Fig.  2-1  that  the  computed  hodograph  gives 
the  wind  distribution  except  in  the  lower  part  of  the  layer  below  the  inversion.  This  defi¬ 
ciency  may  be  explained  as  follows:  (1)  The  topographical  effect  has  not  been  taken  into 
account;  in  this  case  the  terrain  near  Santa  Monica  Airport  tends  to  speed  up  the  flow  near  the 
surface  of  the  earth  as  indicated  by  the  observed  winds  in  the  surface  layer,  however,  the 
average  wind  speed  over  a  flat  plain  in  neighboring  areas  shifts  its  maximum  to  the  middle  of 
the  marine  layer  (Edinger,  1959).  (2)  The  replacement  of  an  inversion  layer  by  a  surface  of 
discontinuity  in  ten^erarure  exaggerates  the  thickness  of  both  the  lower  and  upper  layers. 

The  sea-breeze  velocity,  Vg  ,  computed  from  the  wind  at  500  m  underestimates  its  intensity. 
(3)  K  is  generally  not  constant  with  height  as  assumed. 


It  is  interesting  to  note  that  however  crude  the  model,  it  brings  out  the  characteristic 
features  of  the  stationary  flow  ii,  the  planetary  boundary  layer  with  a  temperature  inversion 
and  a  sea-breeze. 
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5,  Conclusion 

An  analysis  of  the  effects  of  an  inversion  layer  and  a  sea-breeze  on  the  stationary  flow  in 
a  planetary  boundary  layer  shows  that  (1)  an  intense  temperature  inversion  virtually  acts  as  a 
boundary  surface  which  divides  a  planetary  boundary  layer  into  two  layers,  (2)  in  the  lower 
layer  the  flow  is  primarily  thermally  driven  whereas  in  the  upper  layer  the  flow  is  similar  to 
that  in  a  boundary  layer  of  homogeneous  fluid.  The  inversion  layer  moves  in  the  direction  of 
the  resultant  of  the  stresses  acting  on  the  layer  by  the  fluids  above  and  below  it. 
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PART  111 

TURBULENT  TRANSFER  IN  THE  BOUNDARY  LAYER 
OF  A  STRATIFIED  FLUID 
by  S.-K.  Kao 


1.  Introduction 


In  recent  years,  a  series  of  wind  and  temperature  observations  under  various  stability 
cunTdraons-iTi-rbe-surface-tayerdiaj  beernnade-itotabiy-by-Thornthwaite  and  Kaser  (1943),  and- 
by  Deacon  (1949)-  Thornthwaite  and  Kaser’s  observations  brought  out  the  characteristic 
features  of  the  velocity  profile  in  neutral  and  non-neutral  conditions,  whereas  Deacon  found 
thawJje  observed  wind  profile  could  be  represented  by 


with 


u*  k(l-/3) 

H  <  1 
fi  =  I 


for  stable  condition 
for  neutral  condition 


[S  >  \  for  unstable  condition 

✓ 

where  U  is  the  mean  velocity,  u*  is  the  friction  velocity,  k  is  Von  Karman’s  constant,  z  is 
the  height,  Zq  is  the  roughness  parameter,  and  is  a  stability  parameter.  Deacon’s  treat¬ 
ment  has  been  very  useful  in  representing  data  but  it  has  its  limitations.  Firstly,  in  Deacon's 
formula  [i  is  assumed  to  be  independent  of  height.  However,  Davidson  and  Barad  (1956) 
found  from  their  analysis  of  the  micrometeorological  data  obtained  during  Project,  Great  Plains 
that  ft  decreases  with  height  under  stable  conditions  and  increases  with  height  under  unstable 
conditions.  Secondly,  a  more  fundamental  defect  of  Deacon’s  relation  which  has  been  pointed 
out  by  Sheppard  (1958)  is  that  it  can  hardly  arise  from  a  dimensionally  consistent  analysis  of 
the  problem. 

With  regard  to  the  analytic  aspect  of  the  problem,  the  theory  of  turbulent  flow  in  a  non- 
neutrally  stratified  surface  layer  is  notoriously  difficult,  although  the  well-known  logarithmic 
velocity  profile  for  a  homogeneous,  incompressible  fluid  may  easily  be  obtained  either  on  di¬ 
mensional  grounds  or  on  the  assumption  of  a  mixing  lengtii  being  proportional  to  the  distance 
from  the  surface  (Prandtl,  1932).  Rossby  and  Montgomery  (1935)  made  the  first  notable  analysis 
of  the  problem  for  a  stably  stratified  fluid.  Extensive  investigations  of  the  mechanism  of  the 
turbulent  flow  were  made  by  Lettau  (1949)  and  later  by  Businger  (1954).  Applying  the  similar¬ 
ity  principle  to  the  study  of  the  problem  Monin  and  Obukhov  (1954)  greatly  simplified  the  anal¬ 
ysis.  However,  in  their  treatment  the  function  which  would  give  the  effect  of  the  heating  and 
cooling  of  the  surface  on  turbulence  has  not  been  determined.  Various  related  problems  have 
not  yet  been  settled.  Tt  is  the  purpose  of  this  paper  to  analyze  some  characteristics  of  tur¬ 
bulent  transfer,  and  determine  the  velocity  and  temperature  profiles  in  a  stratified  surface 
layer. 
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2.  Formulafion  and  anniysii 

^  ilir  I'!  t»  Mtl  in>Ml  fltllil  •iV'-roil  Idfilllu*  (llliritIMI  lilMJ’lt  V«lltl 

vclm  iiy  v.ui.niDiiy  sniatl  cmnjiKrtd  wltli  tlie  locrtt  vclociiy  ot  sound,  with  temperature  varia¬ 
tions  small  compared  with  the  absolute  temperature,  and  unaffected  by  the  rotation  of  the  earth. 
W'ithout  serious  loss  of  generality,  we  may  suppose  the  mean  flow  to  be  unidirectional  and 
horizontal,  and  describe  the  flow  in  a  Cartesian  coordinate  system  with  the  x  axis  in  the  hori¬ 
zontal  direction  of  the  flow  and  with  the  z  axis  vertically  upward.  In  general  the  vertical 
flux  of  heat  in  the  atmosphere  is  partly  radiative.  For  the  simplification  of  a  theoretical  anal¬ 
ysis  of  the  problem  of  transfer  in  surface  layer  wc  shall  exclude  the  component  of  flux  due  to 
radiative  transfer  and  consider  the  vertical  flux  of  any  entity  to  be  independent  of  height  and 
of  the  horizontal  coordinate. 


by 


In  conditions  of  neutral  stability  it  is  known  that  the  mean  velocity  U  at  height  z  is  given 


(1) 


where  the  subscripts  ^  indicate  the  case  of  neutral  stability,  thus  u*^  is  the  friction  veloc¬ 
ity  in  neutral  condition.  Equation  (1)  implies  that  in  neutral  condition  the  mixing  length  for 
momentum  transfer  1^,^  =  kz. 

Integration  of  (1)  gives  the  well-known  logarithmic  law 

where  the  constant  of  integ^tion  Zq  is  a  length  characterizing  the  roughness  of  the  surface. 
The  validity  of  (1)  and  (2)  for  the  flow  of  air  over  the  land  and  sea  has  been  demonstrated  by 
Sverdrup  (1936),  Thornthwaite  and  Kaser  (1943),  Sheppard  (1947),  Pasquill  (1950),  Rider  (1954), 
Hay  (1955),  and  Deacon,  Sheppard  and  Webb  (1956). 

It  is  obvious  from  (1)  and  (2)  that  in  neutral  condition  the  only  parameter  which  needs  to 
be  considered  is  the  friction  velocity.  The  presence  of  buoyancy  force  and  heat  flux  in  con¬ 
ditions  of  stable  and  unstable.stratifications  suggests  two  parameters:  gT"^  and  Q(pCp)"^ 
where  g  is  the  gravity  acceleration,  T  is  the  mean  temperature,  Q  is  the  vertical  heat  flux, 
p  and  Cp  are  respectively  the  air  density  and  the  specific  heat  of  the  air  at  constant  pressure. 

In  terms  of  these  parameters  we  introduce  the  following  characteristic  temperature  and  length 


,OCp  u* 


kl.  _Q 

T  pc, 


In  a  neutrally  stratified  surface  layer  the  turbulence  is  entirely  due  to  friction.  When  there 
is  a  heat  flux,  turbulence  increases  in  an  unstably  stratified  atmosphere  in  which  convective 
energy  is  released,  whereas  it  decreases  in  a  stably  stratified  atmosphere  in  which  convective 
energy  is  absorbed.  Correspondingly,  the  surface  stress  increases  in  the  former  case  and  de¬ 
creases  in  the  latter.  Friction  velocity  is  therefore  a  functioti  of  heat  flux  and  may  be  ex¬ 
pressed  by 


u*  -u*^  f(l) 
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where  !  ,  c'liiire.r'eri.st'c  length  and  f(hL~*)  is  ..  function  to  be  determined  having  the  prop¬ 
erty  .  /-'■>  (Q  -O',  ■:*  :.  'lu'-e,.  to  u*^,  so  that  f(0)  =  1.  It  is  to  be  noted  that  in  (3)  the 

friction  vt  '.ei.i  y  fc'  neutrai  ..  ability  u*^is  introduced  instead  of  u»  as  introduced  by  Monin  and 
Obukhov  fi '■).■■■  0-  dhis  selection  of  parameter  enables  us  to  determine  later  the  friction  velocity 
in  non-neuodl  conditions  as  a  function  of  heat  flux  and  the  friction  velocity  in  neutral  stability. 

Analogous  to  the  formulation  for  the  flux  of  momentum  and  heat  in  the  laminar  flow,  let  us 
first  introduce  the  concept  of  eddy  viscosity  and  eddy  conductivity  defined  respectively 
by  the  following  momentum-  and  heat-transfer  equation 


K 


m 


dU  ,, 
dz 


u* 


Kh  =  -_Q 

dz  ocp 


(5) 


(6) 


where  6  is  the  mean  potential  temperature.  Now,  on  grounds  of  dimensional  reasoning,  we 
generalize  relationship  (1)  to  conditions  of  non-neutral  stratification 


where  is  a  function  of  the  non-dimensional  variable  z  L  with  the  property  that 
as/L/-®.  Substitution  of  (7)  and  (4)  in  (5)  gives 


which  implies  that  the  momentum  mixing  length 


_  k  z 


(9) 


Similarly,  we  have  for  the  transfer  of  heat 

Ky,  «  kzu*^  _Aki 

^ir) 

dz  kz  f  (1] 


(10) 


(11) 
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Best  Available  Copy 


where  v’]i  L"^) 


is  a  fimccion  to  be  determined. 


Equations  (10)  and  (11)  imply  that  the  mixing 


length  for  heat  transfer 


kz 


(12) 


Since  it  has  been  shown  experimentally  (Richardson,  1920)  and  theoretically  (Batchelor, 
1953)  that  the  Richardson  number  Ri  can  be  used  as  a  measure  of  the  stability  of  the  atmosphere, 
it  is  of  interest  to  observe  from  (7)  and  (11)  that 
g  de  z 

Ri  =  T  ^  =  -1  ~L 

(^2  (13) 

dz  L 


where  y  ~  '/i^l  -  jj^l  .  This  indicates  that  the  Richardson  number  is  gener¬ 

ally  a  function  of  height,  heat  flux  and  eddy  conductivity. 


3.  Acceleration  balance  equation  and  generalized  profiles 

It  has  been  pointed  out  in  the  previous  section  that  atmospheric  turbulence  in  the  surface 
layer  is  caused  by  mechanical  friction  and  thermal  convection.  Total  turbulence  may  therefore 
be  considered  as  the  combination  of  frictional  turbulence  and  convective  turbulence.  Attempts 
at  a  formulation  of  this  relationship  have  been  made  by  Rossby  and  Montgoraer/  (1935)>  Lettau 
(1949)  and  Businger  (1954).  In  their  treatments  these  investigators  implicitly  considered 
frictional  turbulence  as  an  entity  determined  by  the  mean  motion  in  condition  of  neutral  stabil¬ 
ity.  Judging  from  the  nature  of  turbulent  motion  that  the  frictional  turbulence  and  the  total 
turbulence  are  generally  interdependent,  we  shall  consider  the  frictional  turbulence  to  be 
dependent  on  the  mean  structure  of  the  total  turbulence. 

We  shall  formulate  the  frictional,  convective  turbulence  relationship  in  the  form  of  accel¬ 
eration.  It  is  obvious  that  the  total  turbulence  in  the  dimension  of  acceleration  may"be  written 
as  ln,(dU/dz)2.  The  contribution  of  the  frictional  turbulence,  which  gives  rise  to  a  mixing  length 
Im^i  I  acceleration  of  the  total  turbulence  is  ln,^(dU/dz)2.  The  effect  of  buoyancy  force 

on  the  motion  of  an  cd^-  may  be  evaluated  as  follows:  when  an  eddy  has  a  temperature  difference 
T'  with  its  surroundings,  this  eddy  has  a  convective  acceleration  gT'/T  which  on  the  aver¬ 
age  has  a  value  equal  to  -  lh(g/T)U  s/dz).  The  acceleration  balance  equation  for  turbulent 
flow  may  therefore  be  expressed  by 


where  B  is  a  proportionality  factor  as  yet  undetermined.  A  form  similar  to  this  has  been 
studied  by  Elliott  (1957). 

Substitution  of  (7),  (9),  (11)  and  (13)  in  (14)  gives 

-^m  ♦  4-  •  0 

f3  L 
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(15) 


Solvii!"  ior  v'vi;  W'-  (’Iilain 


f^L  4B 


Here  we  select  the  solution  which  satisfies  the  property  -  1  as  /L/-®  ,  which  follows 
from  the  dimensional  reasoning  as  analysed  in  the  previous  section.  Equation  (15)  indicates 
that  i/>j^(zL""^)  has  two  branches,  one  for  upward  heat  fluxes  (Q  >  0)  and  the  other  for  down¬ 
ward  heat  fluxes  (Q  <  0).  In  the  former  0.5  ^  i  1  ,  whereas  in  the  latter  1. 

By  substituting  (15)  in  (7)  and  by  integrating  Zq  to  z  we  obtain 


u*. 


f- 

In  £,  +  (] 

4Bz,H  _  / 

tanh"*/j 

[_4Bz\M  . 

r  tanh-1/1- 

Q>0 

tor 

(2 

z„  ^ 

f3L^  ^ 

f3L  ’  < 

f3L^ 

^  f3L  M 

f3L  4B 

(16a) 

.  In 

z 

for 

Q  =  0  ,  z  > 

(16b) 

f- 

ln£_+  1 

fl-  -  coth-1 

/I-  1153' 

^  +  coth-l/l- 

for  Q  <  0,  z  > 

12 

-0  < 

.  f3L/ 

^  f3L' 

) 

(16c) 


Mathematically  Eqs.  (I6a)  to  (l6c)  are  valid  at  all  heights  with  U  -®  as  z  -®.  In  view 
of  the  fact  that  both  the  sources  of  mechanical  turbulence  and  heat  flux  are  located  at  the  sur¬ 
face  of  the  earth,  turbulence  intensity  must  eventually  decrease  with  height.  Eq.  (16)  applies 
only  to  the  turbulent  surface  layer.  Let  us  assume  that  a  free  stream  velocity,  Ug,  exists  at 
the  height  of  the  thickness  of  the  turbulent  surface  layer,  Hg.  By  substituting  the  free  stream 
velocity  in  (l6b)  and  (16a)  or  (l6c)  respectively  for  conditions  of  stable  or  unstable  stratifi¬ 
cation,  then  by  eliminating  Ug  from  these  equations,  we  may  determine  f  numerically.  In  this 
case  the  height  of  the  thickness  of  the  turbulent  surface  layer  is  selected  as  the  characteristic 
length  h  in  (4).  An  approximate  expression  for  f  explicitly  in  terms  of  heat  flux  will  be  given 
in  the  next  section. 


Rider  (1954)  found  from  the  analysis  of  his  observations  of  wind,  temperature  and  humidity 
that  their  profiles  are  rather  similar  over  the  first  few  meters.  This  implies  that  the  eddy 
viscosity  is  proportional  to  eddy  conductivity.  If  we  assume  that  the  ratio  7 
stant,  Eq.  (11)  may  ^  integrated. 
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Ic  follows  froiii  (13)  :in(.!  (15)  that 


Ri  . 


{‘f-fri 


The  flux  form  of  Richardson  number  (Richardson,  1920;  Deacon,  1955)  Rf,  which  is  equal 
to  the  ratio  of  the  rate  at  which  buoyancy  forces  extract  energy  from  the  turbulence  to  the  rate 
at  which  it  is  supplied  by  the  shear-stress,  may  be  expressed  in  the  following  form 

Rf  =  Ri  ^ 

Substitution  of  (18)  in  the  above  equation  gives 

p, .  ''k 


Richardson’s  arguments  make  It  clear  that  Rf  cannot  exceed  unity  over  any  substantial  depth 
of  fluid,  so  it  must  have  a  maximum  value  (critical  value)  not  greater  than  this.  It  is  readily 
seen  from  the  above  equation  that  Rf  reaches  its  maximum  value  when  4Bz(f3L)  =1.  The 
lower  limiting  value  of  the  flux  Richardson  number  is  therefore 


Rf  =  - 


in  which  B  =  1  is  chosen. 


Townsend  (1958)  has  recently  shown  that  the  flux  Richardson  number  cannot  exceed  0.5, 
which  together  with  (18a)  gives 


The  distribution  of  the  flux  Richardson  number  as  a  function  of  heat  flux  and  height  is  shown 
in  Fig.  3-1.  It  is  shown  in  this  diagram  that  the  larger  a  positive  heat  flux  the  greater  is  the 
rate  of  decrease  of  the  flux  of  Richardson  number  with  height,  whereas  the  smaller  a  negative 
heat  flux  the  greater  is  the  rate  of  increase  of  the  flux  of  Richardson  number  with  height. 
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!n  'r'ig.  ’.-2  the  variiition  of  non-dimensional  mixing  length,  /Zg,  is  exhibited  for  various 
values  of  non-dimensional  heat  fiux,  S  -  ZgkgQfpCpTu^)"!,  the  sloping  straight  line  in  the 
diagram  describes  the  mixing  length  for  the  neutral  case.  The  curves  on  the  left  of  this  line 
represent  the  mixing  lengths  for  stable  cases,  the  curves  on  the  right  of  it  represent  the  mixing 
lengths  for  unstable  cases.  The  slope  of  mixing  length  decreases  with  height  for  a  negative 
heat  flux,  whereas  it  increases  with  height  for  a  positive  heat  flux. 

4.  Small  values  of  ^(f^L)"^ 

For  small  values  of  z(~^  L“^,  'pm  expanded  in 

. forJ^<<l 

L  f3L  f3L  •  f3L  (19) 

By  substituting  (19)  in  (7)  and  (11),  then  by  integrating  from  z  to  Zg,  we  obtain  respectively, 


U(z)  =  rin^ 


z  _  B 


Zo  f3L 


(z  -  Zn 


0{7.)  =  s(Zg)  +_il  [in  i_-  ®  (z  -  Zg)  T 
kfyl-  z„  f3L  J 


This  gives  the  (logarithmic  +  linear)-law,  the  same  result  obtained  by  Monin  and  Obukhov 
(1954)  through  the  expansion  of  an  undetermined  function  with  an  empirical  value  of 
B  =  0.6  ±0.06.  Monin  and  Obukhov’s  findings  are  therefore  the  first  approximation  of  the  veloc¬ 
ity  and  temperature  profile  (16)  and  (17). 


To  determine  f,  let  us  substitute  the  free  stream  velocity  U3,  at  height  Hg,  in  (20), 


which  becomes  for  neutral  stability  (/L/-®  and  f  =  1), 

Ug  =  ^  In  (23) 

Here  we  assume  that  the  free  stream  velocity  is  unaffected  by  the  heat  flux  in  the  surface 
boundary  layer.  Elimination  of  Us  from  (22)  and  (23)  results  in 

f3  -  f2  -  1(3.'..^°^  .  =0  (24) 

L  In  (Hs/zg)  ^  ^ 


Although  it  can  readily' be  shown  that  the  real  solution  of  (24)  is 
f  =  1  +  i  i  i+JZf  ri+  (1+  ^  3  i+I?  f  r  1. 
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where  ^  B(Hs  -  ^-q)  [L'n(Ifs/^o)^  ~^>  expression  for  f  simpler  than  (25)  may  be  obtained. 
This  will  simplify  the  computation.  A  careful  examination  of  (16)  indicates  that  at  a  certain 
height  between  and  Hg  the  velocity  shear  in  non-neutral  stratification  is  equal  to  that  in 
neutral  stability.  Let  this  height  be  hg;  we  have  by  definition 


dU 


dz 


z  =  ho 


z=ho 


(26a) 


Substitution  of  (1)  and  (7)  in  the  above  equation  with  the  use  of  (15)  gives 


Solving  for  f  we  obtain  ^ 


f2  _  £  _  =  0 

L 


(  =  'A  [u(i* 


Here  we  select  the  solution  which  satisfies  the  condition  that  f  =  1  as/L/-aj. 
It  can  be  shown  by  substituting  (1)  and  (7)  in  (26a)  with  the  use  of  (19)  that 

f3  .  £2  _  Bho  ,  0 

L 

Comparison  of  the  above  equation  with  (24)  gives  the  value  of  hg 

He  -  Z„ 


ho  = 


‘s 


By  substituting  (28)  in  (27)  we  obtain 


f  -  Hri  +  (1+  ~  )^-j  -  1  + 

L  Lln(H3/Zg)  i 


B(Hs  -  Zg) 
Lin  (Hg/Zg) 


It  follows  from  (4)  that 


u*  =  u*  11+  1 

I  InfH  /z  'll 

'  *•  ""  '-*S'  O'  ' 


(26b) 


(27) 


(28) 


(29) 


This  indicates  that  an  upward  heat  flux  (Q  >  0)  increases  the  friction  velocity,  whereas  a 
downward  heat  flux  (Q  0)  decreases  the  friction  velocity.  This  agrees  with  what  we  have 
anticipated  from  the  energy  consideration  in  section  2. 

Although  for  small  values  of  z  L~^  velocity  profiles  may  be  estimated  by  (22),  a  more 
precise  distribution  of  non-dimensional  velocity  U/u^^  (16)  is  shown  in  Fig.  3-3  in  which 
Hg/Zg  “  3.10^  is  chosen.  The  sloping  straight  line  in  the  diagram  describes  the  logarithmic 
profile  for  the  neutral  case.  The  curves  on  the  left  of  this  line  represent  the  velocity  profiles 
for  stable  cases;  the  curves  on  the  right  of  it  represent  the  velocity  profiles  for  unstable  cases. 
At  the  roughness  height  the  velocity  shear  in  stable  conditions  is  smaller  than  that  in  the 
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nfMitrn!  condition,  whereas  the  latter  is  smaller  than  that  in  unstable  conditions.  A  character¬ 
istic  feature  of  these  p.'ofilcs  is  that  the  velocity  profiles  for  the  non-neutral  cases  converge 
toward  cite  profile  for  the  neutral  case  in  the  higher  level.  This  agrees  with  the  classical  pro¬ 
files  observed  by  Thorathwaite  and  Kaser  (1943). 


It  is  of  practical  interest  to  determine  the  velocity  profile  and  heat  flux  from  wind  observa¬ 
tions  alone.  For  example,  if  we  have  wind  observations  at  three  levels,  zj,  Z2  and  zj,  we  can 
determine  the  velocity  profile  and  heat  flux  as  follows:  with  the  use  of  (20)  it  can  be  shown 
that 


u*  =  k  (^2  ~  ~  ^^3  ~  ^2^  ^^2  ~ 

(Z3  -  Z2)  In(z2/zj)--(z2  -  zj)ln(z3/z2) 


(31) 


f^L  =  B  (^2  ~  ^1)  (^3  ~  ^2)  ~  (^3  ~  ^2^  (^2  ~ 

(U2  -  Ij2)ln(z3/Z2)-(U3-  U2)ln(22/22) 


(32) 


By  substituting  (31)  and  (32)  in  (20)  we  can  determine  the  wind  profile  in  the  turbulent  surface 
layer.  The  hert  flux  can  be  determined  from  (31),  (32)  and  (3  ), 


pCpk^T 

SB 


[(U2-Ui^1n(z3/22)-(U3"U2)ln(z2/2j)],[(U2-Ui)(z3-22)-(U3-U2)(z2-Zi)]2 


[  (Z3-Z2)  (^2/21)  -  (22-21)  In  (Z3/Z2)  ]  3 


If  both  wind  and  temperature  observations  are  available,  velocity  and  temperature  profiles 
may  be  determined  from  observations  at  two  levels.  It  can  be  shows  with  the  use  of  (20)  and 
(21)  that  ;  .  ' 

,  (U2-Ui)2  -  b|-  (^1-  «2)(22-z.i) 

(U2  -  up  In  (22/21) 


Q  =  pCp  u» 


^1  ~  ^2 

U2-U1 


(35) 


Here  7=1  is  assumed.  By  substituting  (29)  in  (34)  the  frict  ion  Vfflocity  for  neutral  condition, 
u,^^,  may  be  determined. 
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X  Conclusion 


An  analysis  is  made  of  some  characteristics  of  the  steady  turbulent  transfer  in  the  bound¬ 
ary  layer  of  a  stratified  fluid.  The  effects  of  the  heat  flux  on  the  variation  of  the  mixing  length 
and  the  flux  of  Richardson  number  are  determined.  The  velocity  and  temperature  profiles  are 
derived.  It  is  found  that  velocity  profiles  for  non-neutral  cases  converge  in  the  higher  level 
towards  the  profile  for  the  neutral  condition,  a  characteristic  which  agrees  with  velocity  pro¬ 
file  observed  by  Thornthwaite  and  Kaser  (1943).  It  is  shown  that  Monin  and  Obukhov’s  findings 
become  the  first  approximation  of  equations  (l6a)  and  (I6c).  For  a  constant  free  stream  veloc¬ 
ity  an  upward  heat  flux  increases  the  friction  velocity,  whereas  a  downward  heat  flux  decreases 
the  fricion  velocity.  The  lower  limiting  value  of  the  flux  Richardson  number  is  found  to  be 
H3.5,  which,  together  with  the  upper  limiting  value,  0.5,  obtained  by  Townsend  (1958),  give.s 
the  range  of  the  flux  Richardson  number. 
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